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We consider the spin and pseudospin (charge) response functions of the ex- 
actly soluble Anderson atom model. We demonstrate, in particular, that a 
deviation from the magnetic Curie-law behaviour, appropriate for a free spin 
one-half, increases with increasing asymmetry and temperature. In general, 
oscillator strength is transferred from the spin degrees of freedom to the pseu- 
dospin modes. We also consider the negative-U Anderson atom and demon- 
strate that the pseudospin modes are the relevant low-energy excitations in 
this case. Especially, the roles of the spin and charge excitations are inter- 
changed upon reversal of the intrasite Coulomb repulsion, U. 
PA CS numbers: 71.10.-w, 71.27. -ha, 73.21.La, 75.75. -ha 



1. INTRODUCTION 

A single quantum dot behaves like an artificial atom when electronic 
confinement in the dot approaches atomic size.^"^ In such structures, there 
emerge the characteristic features of an atomic impurity: the quantization 
of charge and energy. It is appropriate to describe these artificial atoms by 
means of the Anderson impurity model. "^'^ 

The Anderson model was first proposed to describe magnetic impuri- 
ties in a metal. ^ In the Anderson model, the nonmagnetic-magnetic tran- 
sition of the local d-state was first described within the Hartree-Fock (HF) 
approximation,^ which amounts to truncating the model Hamiltonian into 
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a bilinear form in the fermion operators. This approximation leads to an 

abrupt phase transition, whereas the actual change must be a gradual one for 
a finite system. Numerous ingenious approaches to describe the complicated 
many-body problem associated with the emergence of the interacting cor- 
related many-electron problem have been introduced, such as Green's func- 
tion methods,^ functional-integral techniques,^ numerical renormalization 
schemes,^' ^'^ Bethe-Ansatz approaches'^ and noncrossing approximations.'^ 
Perturbation theories have been utilized starting both from the small-^7 (HF) 
limit'^"'^ and from the atomic (V^ = 0) limit. '^''^ Furthermore, interpola- 
tion schemes between the HF and atomic limits have been developed. '^"^' 
However, the properties of the atomic limit, which again is exactly soluble, 
appears not to have been thoroughly discussed in the literature.^^ 

In this paper we want to discuss in detail the coupling of the correlations 
which takes place in the magnetic (zero-temperature susceptibility is diver- 
gent) Anderson atom. Relevant information on the correlations among the 
electrons may be obtained by considering the response of a system to an ex- 
ternal perturbation. We derive the static response functions of the Anderson 
atom. Our derivation demonstrates that it is in general important to dis- 
cuss not only the spin susceptibilities (x) but also the charge (pseudospin^^) 
susceptibilities (^); the superscript zero in x° and denotes static response. 

We show that the Anderson impurity atom best follows the magnetic 
Curie law = (42^)"', appropriate for a quantum-mechanical spin one- 
half, in the symmetric situation at low temperatures (T <C U). At high 
temperatures, on the other hand, the spin and charge degrees of freedom 
become equally important and the respective susceptibilities approach the 
common limiting value T^^ = Typ ^1/8 for T ^ oo in zero field. 

In particular, our equations may be applied to negative-[/ situa- 
tions. The negative-C/ Anderson model was first proposed to describe 
the electronic structure of amorphous semiconductors.^^ Since then, 
negative-?/ behaviour has been observed, e.g., in the context of high- 
tcmpcrature superconductors,^^ hcavy-fcrmion systems^^ and interstitial de- 
fects in semiconductors.^'^' It has been observed that in a quantum dot, 
the second electron in the dot may be more strongly bound than the first one 
under certain circumstances. ^'^"^^ Essentially two different mechanisms have 
been proposed to supply the effective net attraction between the electrons to 
cause a quantum dot negative-C/ properties. ^^'^^ However, the microscopic 
origin of the phenomena is still unclear. 

The reversal of U changes the roles of spin and charge. For the symmet- 
ric level configuration and negative f/, the charge response functions behave 
like the spin susceptibilities for positive U . Hence, for C/ < 0, the spin de- 
grees of freedom are frozen out and the charge degrees of freedom are the 
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dominating low-energy excitations. In what follows, we intend to elucidate 
the dual roles played by the magnetic field in connection of spin dynamics 
and that of asymmetry in the context of charge dynamics. 

For high magnetic fields, we find new peak structures in the longitudinal 
spin and charge susceptibilities and an associated threshold behaviour in the 
corresponding transversal response functions. These features are associated 
with level crossings occurring between states belonging to different Fock 
spaces in special configurations displaying high symmetry. 

The rest of this paper is organized as follows. The Anderson Hamilto- 
nian is presented in Section 2. In Section 3, we take a close look into the spin 
and charge algebras. Section 4 is devoted to derivation of analytic expres- 
sions for the spin and charge susceptibilities. In Section 5, we present our 
numerical results. In Section 6 we give discussion and conclusions. In Ap- 
pendix A we collect the central properties of double-time Green's functions 
that are utilized throughout the present paper. 

2. MODEL HAMILTONIAN 

The Anderson Hamiltonian for magnetic impurities in metals^ 

W = X] ^karika + {^kcl^^a + V^dlcka) + ^ E^n^ 

k,(T k,a f 

+Un^n^ (1) 

describes the transition of the local d-electron orbital from a nonmagnetic 
resonant virtual bound state (F/C/ » 1) to a magnetic atom (F/C/ <C 1). 
Here U is the intra-atomic Coulomb repulsion energy and F = iTN{Qi){\Vk\^) 
(with A^(0) the density of conduction electron states at the Fermi level, 14 
the d-level hybridization matrix element and with () denoting an average 
over the Fermi surface) is a measure of the admixture of the local state with 
energy = E — aB (here cr = ±1/2 and B is the external magnetic field). 
Furthermore, c\.^, c^.^-, dj. and d^j arc the creation and annihilation operators 
for electrons in the conduction band and in the impurity state, respectively, 
and the corresponding occupation-number operators are n^^ = c^^Cka and 
ricr = d]jd(j. The conduction-electron dispersion relation is denoted with Sk^. 
In the nonmagnetic (d-spin susceptibility is finite at T = 0) [/ = limit, 
the Hamiltonian (1) is of bilinear form and hence the static and dynamic 
properties of the impurity spin may be obtained exactly in closed form."^^ In 
this paper we discuss the properties of the Anderson model in the atomic 
{Vk = 0) limit 

7ia = T.aE„n„ + Un^n^, (2) 
for which only the last two terms in the Hamiltonian (1) remain. 
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Here we form the spin and charge operators S and Q. We also sum- 
marize some of their properties. We define the operator-valued spinor 

ip = {dif,di)'^, and introduce spin-operator components as Sa = ip^cTaip 
for a G {x, y, z} and = {Sx ± Sy) /a/2- Here aa are the usual Pauli spin 
matrices. Thus we obtain 



S+ = —d\d^ (3) 



1 

S- = ^4^T (4) 
= sl = S^y = Un^+n^-2n^ni) = \s\ (6) 



With the help of the spinor (f) = (^d^ , c/j j , the charge-operator components 
may be expressed analogously in the form Qq, = (f>^aa(t>- We find 



= (7) 

Q' = -^d^d^ (8) 

Qz = lin^+H-l) (9) 

Ql = Ql = Ql = ]{l-n^-n^ + 2n^n^) = \Q\ (10) 



From the operators S and Q, we can easily form a spin-1/2 algebra C which 
obeys the canonical spin commutation relations [Ci,Cj]_ = if^jkCk with 

= c{c + 1) for c = 1/2, as follows. We note that S and Q obey by 
construction the canonical commutation relations. Moreover, for the atomic 
Anderson Hamiltonian we find that Sz and Qz and thus Sl and Ql (therefore, 
also S*^ and Q'^) are constants of motion. However, {S"^) = | is not obeyed, 
unless (Q^) = 0, and vice versa. Furthermore, one easily finds that Qj\ = 
for all i, J. Consequently, a spin-1/2 algebra C can be formed as C = S+Q. 
One easily sees that = S"^ -I- = |, as required for a spin of fixed 
magnitude 1/2. This clearly demonstrates the importance of considering 
both spin and charge degrees of freedom. 
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4. EXPRESSIONS FOR THE SPIN AND CHARGE 
SUSCEPTIBILITIES 



38 



By utilizing the Hubbard operators 



A^^ = {l-n^^)da (11) 

A2a = n-ada (12) 

(for a recent discussion of supersymmetric Hubbard operators, see^^) or the 
spin-flip and charge-transfer operators 

= dl^d^ = V2S^ (13) 
Ai^ = d-ada = V2Q^ (14) 

the atomic Hamiltonian (2) may be represented as a bilinear form in two 
alternative pictures 

{aIa,, + aIA2.) + UAl^A2a (15) 

= S^E^(4^A3^+4^A4^) +C/4^^4a- (16) 

Action of the operators Ai^, A^a, A^a and A^^ and their hermitean adjoints 
on the states |0), |t), |i) and |ti) is shown in Table I. These operators are 
also eigenoperators of the Hamiltonian 

[Aka, Ha] = aka{E, U, B)Aka, (17) 

where aka{E, U, B) is a scalar- valued function of the model parameters. Con- 
sequently, it is easy to obtain the following anticommutator (+) functions 
directly from the equations of motion (46) and (47) (no coupling to higher- 
order Green's functions) 

((^i.;4.))^ = ^^^f^ (18) 

{{A^.;Al))t = i^^ + ^-^-J^'^^-'r) (20) 

The corresponding commutator (— ) functions are found by replacing the 
expectation values in Eqs. (18)-(21), respectively, with ((1 — n_o-)(l — 2no-)), 
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Table I 

Operators Aiu, A2a, ^So- and and their Hermitean adjoints describe all 
the possible transitions between the states jO), |t), ||) and |||). Notation 
means that one obtains the final state |/) through the operation of an entry 
in the table on the initial state 



initial 



final 



I/) 



|0) 

IT) 
li) 
ITi) 



|o) IT) li) ITi) 



1 

4 



T 



i 

^3T 



A 



n 



Al 



31 



A^ 

a\-^ a. 



Ak 



21 ^2T 



A41 

M 
A21 
i 



{n-.fj — 2n(jn^(j)^ {n^fj — n„) and {l—na—n^fj) . Note that for the d electron we 
have d(j = Aifj+A2a and for the d-electron propagator we find: {{da ', d]^))~^ = 
{{A,a ; Al))t + {{A2a ; 4,))+ = (l-n_.)/ (z - Ea) + {n.a)/ {z - - U) , 
with poles at the single-particle eigenenergies of the atomic Hamiltonian (2), 
z = E(j and z = E^ + U . 

Considering the time correlation functions (50) and (51) for the Green's 
functions (18)-(21), one finds after some algebra 



{ria) 



{naU. 



where 



{naU-a) = 



f{E_a + uy 

- f{Ea) - f{E_a) 



1 



f{Eo 



f{E^ 



(22) 



(23) 



fi2E + U) f{Ea + U) f{E_a + U) 
is the correlated double occupancy and / is the Fermi function. 

Now we turn our attention to the longitudinal and transversal spin and 
charge response functions of the Anderson atom. They may be defined, for 
complex frequencies, as 



X±(^) = -((5+;5-)),- = -i((^3i;4i))." 

^±{z) ^ -((Q+;Q-))J = -^((4^;A4t))J 

X\\{z) = -{{[s,-{s,))-[s,-{s,))))- 

^II(Z) ^ -(((q,-(Q,));(q,-(Q,))))-. 



(24) 

(25) 

(26) 
(27) 
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The expectation values {S"^) and (Q^) vanish for Ha and are thus not needed 
in the defining equations (24) and (25). However, (Sz) and (Qz) in Eqs. 
(26) and (27) are generally nonzero. In connection with the 0(3)-symmetric 
Anderson model and the two-channel Kondo model, the charge operators are 
called isospin operators^*^ and the parallel charge susceptibility is defined as 
a Qz — Qz response as above but without extracting the mean. 

The transversal spin susceptibility ((— ) function) is readily found from 
Eq. (20) (changing the expectation value as described in the text after (20)) 

xx(.) = (28) 

Consequently, the familiar static zero- frequency limit of x±(z) is found with 
(49) 

xi=r^^^=%^. (29) 

Similarly, for the transversal charge susceptibility one obtains 



and 



/ di^ ii_{^) {Qz) 
^^=/ -—r = -2EV(J- ^^^^ 



oo 



Here 2E + U measures the asymmetry of the level configuration with respect 
to the Fermi level, vanishing in the symmetric {E = —U/2) situation. Thus 
asymmetry behaves for the charge degrees of freedom as the magnetic field 
for the spin degrees of freedom (compare Eqs. (29) and (31)). 

For the parallel spin response in Eq. (26), one cannot directly calculate 
the commutator (— ) function with the equations of motion (46) and (47) 
since [Sz — {Sz),Sz — (Sz)] = 0, and also [Sz — {Sz),Ha] = 0. However the 
corresponding anticommutator function is easy to find since {Sz — {Sz), Sz — 
{Sz)} = 2{{Sl) - {Sz)"^), after which the static parallel spin susceptibility 
Xy is obtained with the help of the fluctuation-dissipation theorem (54). The 
result is 

X8 = i^Slz^)!. (32) 
For the parallel charge response, a similar calculation yields 

^0 _ {Qz) ~ {Qzf' ^23^ 

The parallel susceptibilities may, furthermore, be calculated from 

y° - ^ (34) 
^11 - dB ' ^^^> 
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and 



(35) 



d{2E + U) 

where the partial derivation is performed such that the intra-atomic Coulomb 
repulsion U is held constant. 

The expectation values {Sz), {Sf), {Qz) and (Q^) in the above results for 
the response functions are expressible through {ria) and {n^n-u) as shown 
in Section 3. Furthermore, since we know {ua) and {naU-a) from Eqs. (22) 
and (23), respectively, it is now straightforward to find the susceptibilities. 
In terms of the dimensionless parameters 



X 

S 



T/U, y = E/U, b = B/T, 



the final results are 



(36) 



t4 = 



sinhl 

1 



b 
2 

- coshl 
2 \2 



{hS) 

b 



|e~-^/^ + ^cosh^ 



[2{l + 2y)S] 



(37) 

(38) 
(39) 

(40) 



In what follows, we omit the superscript zero from y and y since 
in this paper we only consider the static response functions. 

In addition, we mention the following properties obeyed in the zero-field 
and symmetric limits for the static responses 



X|| = lims-^o X± (41) 
i;i™tr/2^ll =l™B-.-t//2C±- (42) 

The equalities (41) and (42) again demonstrate that the asymmetry is anal- 
ogous to an external magnetic field when one considers the charge response 
functions instead of the spin susceptibilities. All the response functions are 
symmetric with respect to reversal of the external field {B —B). What 
is more interesting, however, is that the susceptibilities are also symmetric 
with respect to reversal of asymmetry (a = 2E + U ^ —a). 
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5. RESULTS 

Here we consider the magnetic-field and level-asymmetry dependencies 
of the spin and charge response functions, X\\,±{T) and ^||,_l(T'). We keep 
U constant and vary E (or y), such that there is a one-to-one relationship 
between y = E/U and asymmetry a = 2E + U = U{2y + 1). Consequently, 
the curves are seen to be drawn for varying level asymmetries. Figures 1- 
7 are for ?7 > while Figs. 8-12 describe the negative-C/ situation. The 
dashed curves always denote y < 0, while the solid ones are for y > 0. 
Furthermore, we have labelled the curves such that the field parameter b = 
B/T is shown inside the parentheses. When the field parameter is omitted, 
6 = is implied. Owing to the above-mentioned symmetry of the response 
functions with respect to reversal of asymmetry, it is convenient to consider 
positive asymmetries only for J7 > and negative asymmetries for U < 0. 
In particular, this means that y > —1/2 below. 

The number of illustrations that follow is large. However, we explore 
the four response functions to illustrate their interrelationships and symme- 
tries upon reversal of U with varying level asymmetries and external field 
strengths. 

5.1. Positive U 

5.1.1. Zero External Field 

Figure 1 shows the spin susceptibility in zero external field (x_l = X||) for 
various asymmetries. This is also given in reference.^ We observe that the 
Curie law for a quantum-mechanical spin one half, Tx = 1/4, is best obeyed 
in the low-temperature limit for the symmetric situation (curve "a" ). How- 
ever, for y < 0, the spin susceptibility always finally rises to the level {AT)~^ 
for low enough temperatures. This happens simultaneously with the depres- 
sion of the parallel charge susceptibility, see Fig. 2. Also the perpendicular 
charge susceptibility vanishes at low temperatures for all configurations, see 
Fig. 7a. It is thus legitimate to state that the spin degrees of freedom are 
the relevant low-energy (or strong-?/) excitations for [/ > 0. Furthermore, 
we find for zero field {{Sz) = 0) 



Tx+Ti\\ + {Qzf = Tx+Ti\\ + 




Therefore, as the charge susceptibilities vanish for low temperatures, we 
obtain Tx = 1/4, as already stated. For high temperatures (T S> i7), on the 
other hand, the oscillator strength becomes evenly distributed among the 
spin and charge degrees of freedom and the spin and charge susceptibilities 
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B=0; U>0 
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Fig. 1. Curie parameter Tx in zero external field (x \\ = X±) foi" positive U as 
a function of level asymmetry. Here \y\ G {0} IJ | p G {1,4, 7, 10, 13, 16}}, 
such that for the curves "a" ^"f ("g" ^"1" ) , y is negative (positive). The 
curve labeled "m" denotes y = 0. 

approach the common value ^"i-^T/u ^oo'^X = l™T/;7^oo^C = 1/8 in zero 
field for all values of y, c.f. Figs. 1, 2 and 7a. 

The y = configuration is a special case. From the general result for the 
spin-state occupations in Eq. (22), it is easy to sec that here the occupation 
numbers satisfy: (n|) = (n|) = 1/3, which yields {Qz}'^ = 1/36, see Eq. 
(10). Consequently, one finds that the low-temperature limits in zero field 
for y = 0: Tx = 1/6, T^y = 1/18 and T^± are consistent with the 
result in Eq. (43). 

5.1.2. Finite Field 

Figures 3a and 4a show the field dependencies of the longitudinal and 
perpendicular spin susceptibilities, respectively, for weak fields. One ob- 
serves that the longitudinal component is more strongly affected by the 
external magnetic field. Furthermore, in high fields - Figs. 3b and 4b - 
the longitudinal spin response displays peaks of invariant height for y > 0, 
whereas the transversal function exhibits threshold behaviour. These peaks 
and thresholds correspond to the situation where one of the localized energy 
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B=0; U>0 




1e-6 1e-4 0.01 1 100 

T/U 



Fig. 2. Curie parameter for the charge susceptibihty T^n with U > for van- 
ishing magnetic field as a function of asymmetry. Here |y| G {0} [j{-^\p & 

{1, 5, 10, 16}}, such that for the curves "a" ^"d" ("e" ^"h" ) , y is negative 
(positive). The curve labeled "i" denotes y = 0. 
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(a) U>0 (b) U>0 
, 1 7/32 1 ^ 




Fig. 3. Magnetic-field dependence of the longitudinal magnetic susceptibility 
X\\{B,T) in (a) moderate and (b) extreme fields. The curves "a" and "b" 
are for \y\ = 1/2 while "c" and "d" denote \y\ = 1/(2^^). The curves for 
which y = are labeled "e" . Magnetic field is measured in terms oib = B/T 
which is indicated in the parentheses in this figure. 




levels crosses the Fermi level as shown in Fig. 5a for y = 1/2. The height of 
the peaks in Tx\\ is 1/16, irrespective of the absolute value of y. 

The parallel and perpendicular charge response functions for weak and 
strong fields are considered in Figs. 6 and 7. Also for the charge suscepti- 
bility, the parallel component is peaked at high fields for y > at the level 
1/16. This is shown in Fig. 6b for y = 1/2 in which case the peak again oc- 
curs in the level configuration of Fig. 5a. Simultaneously, the perpendicular 
charge susceptibility shows threshold dependence, see Fig. 7b. In Fig. 7a, 
for T^± the curve for y = would lie between the |y| = 1/16 curves shown 
in the figure. 
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Fig. 5. Special values of large applied magnetic fields yield level configura- 
tions with particular symmetry. These level schemes correspond to distinct 
features in the magnetic and charge response functions (a: c.f., Figs. 3b, 4b, 
6b, 7b; b: c.f.. Figs, lib, 12b; c: c.f. , Figs. 9b, 10b, lib, 12b). 
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(a) U > 



(b) U > 









d(4) d 

e r\/' 












/ '■' \ ' ' ^ 
8(4)/ / 













b(256) b(32) 



Fig. 6. Field dependence of the parallel charge susceptibility in (a) 
moderate and (b) extreme fields. The curves "a" and "b" are for \y\ = 1/2 
while "c" and "d" denote \y\ = 1/(2^*^). The curves for which y = are 
labeled "e" . The zero-field situation was already considered in Fig. 2. 



(a) U>0 (b) U>0 
^ 1 1/16| ' 




Fig. 7. Field- and asymmetry dependence of the perpendicular charge sus- 
ceptibility ^± foiU>0 in (a) moderate and (b) extreme fields. The curves 
"a" and "b" are for \y\ = 1/2 while "c" and "d" denote \y\ = 1/16. 
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B=0; U<0 




1e-6 1e-4 0.01 1 100 

-T/U 



Fig. 8. Asymmetry dependence of the longitudinal zero-field charge sus- 
ceptibility, ^\\(B = 0,r), for U < 0. Here \y\ G {^}[J{§^\P e 
{16,13,10,7,4,1}}, such that for the curves "a" ^"g" ("h" ^"n" ), y is 
negative (positive). 

5.2. Negative U 

The charge susceptibility in zero field and for weak and strong magnetic 
fields is illustrated in Figs. 8 and 9 for the parallel component and in Fig. 10 
for the perpendicular part. The low-temperature charge Curie law = 1/4 
is here found to be obeyed only in the symmetric situation (y = —1/2, 
C± = irrespective of the magnetic field. Furthermore, an infinitesimal 
asymmetry is sufficient to depress the low-temperature limit as shown in 
Fig. 8a, where the curve "b" corresponds to y ~ —0.499992. Again, in high 
fields a peak of height 1/16 is formed in the parallel charge response T(^jj for 
y > 0, while the perpendicular component displays a threshold at the same 
point. However, here for y = 1/2 this behaviour corresponds to the level 
configuration shown in Fig. 5c, where B = —3U. 

For U < the spin susceptibilities behave in a similar way as the charge 
response functions for U > 0, freezing out for low enough temperatures, as 
shown in Figs. 11 and 12. However, the y = situation is not a special case 
here. Again we observe that the parallel spin susceptibility is more sensitive 
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(a) U < 



(b) U < 
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c,c(4) 














y J" 






Fig. 9. Field dependence of ^y, for [7 < in (a) moderate and (b) extreme 
fields. The curves "a" and "b" are for \y\ = 1/2, "c" and "d" denote 
\y\ = 127/256 and the curves for which |y| = 1/4 are labeled "e" and "f . 




Fig. 10. Field dependence of ioi U < in (a) moderate and (b) strong 
fields. 
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(a) U < 



(b) U < 





1/16- , b(1024) b(256)a(32)b{32) ^ 

a(1024} 



Fig. 11. (a) For U < 0, the longitudinal spin susceptibility x\\ is strongly 
suppressed for T < U/10. (b) A sharp peak is formed for extreme fields 
for both positive and negative values of y. Curves "a" and "b" are for 
|y| = 1/2, corresponding to level configurations in Figs. 5b and 5c, while the 
other curves in (a) represent \y\ = 1/2^^. 



to the external field than the perpendicular component. Furthermore, we 
find peaks in Tx\\ and thresholds in Tx± for high fields, but now both at 
y > and y < 0. The y = 1/2 peaks and thresholds in Figs, lib and 12b 
correspond to the level configuration of Fig. 5c (B = —3U), while those for 
y = —1/2 correspond to the situation in Fig. 5b, where B = —U. 

5.3. Contrasting Spin and Charge 

As pointed out above, the reversal of U exchanges the mutual roles 
of the spin and charge degrees of freedom. We find, in particular, that at 
zero field for the symmetric situation {y = —1/2): ^{U < 0) = x{U > 0) 
and ^(?7 > 0) = x{U < 0). This is illustrated in Fig. 13 for the parallel 
susceptibilities. Furthermore, we find that the spin and charge response 
functions are symmetric with respect to the level asymmetry a = 2E + U, 
such that x(a) = x(— a) and ^(a) = ^(— a) are obeyed. 



6. DISCUSSION 

The Anderson impurity model has proven capable of describing a re- 
markable variety of different physical systems in the field of strongly corre- 
lated electrons. Presently, an important trend in condensed-matter physics 
and nanoelectronics is one where artificial man-made objects are studied, 
rather than real atoms or molecules. The progress in lithography techniques 
has made it possible, e.g., to fabricate quantum dots with properties simi- 
lar to those of real atoms. Our study has relevance for such systems where 
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(a) U < 



(b) U < 





Fig. 12. (a) For U < 0, the transversal spin susceptibility x± is strongly 
suppressed at low temperatures, (b) Extreme field values produce a thresh- 
old for both y > and y < corresponding to the peaks in Fig. 11. The 
curves have been drawn for the same asymmetries as in Fig. 11. 



isolated localized energy levels occur. It would be interesting to investigate 
whether indications for negative-?/ behaviour previously found, e.g., in the 
context of high-temperature superconductivity and heavy-fermion systems 
could also be experimentally realized in quantum dots^^ or carbon nanotubes 
with magnetic impurities. 

We have studied the spin and charge susceptibilities of an artificial An- 
derson atom for arbitrary values of the model parameters. General results for 
all the four relevant response functions have here been presented to the best 
of our knowledge for the first time and an extensive survey of the properties 
of these susceptibilities was carried out. It was pointed out, in particular, 
how the level asymmetry behaves for the charge degrees of freedom as the 
magnetic field for the spin, and that the reversal of U changes the mutual 
roles of spin and charge. In particular, for low temperatures and with pos- 
itive U, the Curie law for a free spin-1/2 is followed by the spin response, 
whereas the charge excitations become suppressed. On the other hand, for 
negative U at low temperatures, the charge susceptibility follows the Curie 
law in the symmetric case, whereas the spin responses vanish. At high tem- 
peratures and for increasing asymmetry, the transfer of oscillator strength 
from the spin degrees of freedom to the charge modes becomes increasingly 
important. In the T ^ oo limit, thermal fluctuations average over any de- 
tails in the energy-level structure. Consequently, the spin and charge modes 
have equal oscillator strengths and they become equally relevant. 

It is noted that for increasing magnetic-field strength, the longitudinal 
susceptibilities are more strongly suppressed than the perpendicular com- 
ponents. Furthermore, at extremely high fields, the longitudinal responses 
show peaks of invariant height while the perpendicular responses have a 
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Fig. 13. Contrasting spin (x||) (solid line) and charge (^y) (dashed line) 
response functions for [/ > (a and b) and for U < (c and d). Note, 
in particular, the spin-charge duality: exchange in the roles of x and ^ for 
U —U (between a and d). For further discussion, see the main body of 
text. 
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threshold. These features are associated with particular level crossings, as 
discussed in the text. 

Our thorough discussion of the spin and charge (pseudospin) suscepti- 
bilities in the atomic limit of the Anderson model has revealed several new 
and interesting details in these response functions. These new features are 
expected to be particularly relevant in connection with the behaviour of 
quantum dots since in this case it is possible to create very high effective 
magnetic fields far beyond those encountered in atomic physics. Therefore, 
quantum dots serve as interesting laboratory models for Anderson model 
physics in limits which have not been explored before. Our results serve 
to emphasize and increase understanding of the relationships governing the 
spin-charge duality in the atomic limit of the Anderson model. The atomic 
model can be taken as the starting point of perturbation expansions in the 
Schrieffer- Wolff limit where the Anderson model can be related to the Kondo 
model for magnetic impurities in metals. Unified definitions of spin and 
charge susceptibilities have been given and they will be utilized in future 
works. 



In order to describe the response of an operator A due to a perturbation 

coupled to operator it is natural to investigate the generalized admittance 
function,^ defined in Zubarev's^^ notation as 



which is a complex function of the frequency variable z and has been an- 
alytically continued for complex arguments. The upper and lower signs 
designate Im(z) > and Im{z) < 0, respectively, and correspond to the 
retarded (analytic in the upper half of the complex plane) and advanced 
(analytic in the lower halfplane) functions. The superscript (+) refers to an 
anticommutator ([A, = {A,B}, correlation) function, while a commu- 

tator ([^, = -B], response) function is meant with the superscript 
( — ). Do not confuse (it) (signs in parentheses) marking the commutator 
(— ) and anticommutator (+) functions with ib or =F (without parentheses) 
denoting the retarded and advanced functions in (44). 

The time evolution of operators is ruled in the Heisenberg representation 
by the operator equation 



A. RESPONSE FUNCTIONS 



{{A-B))f^=Ti dte-*0(±i)([^(i),B(O)](±)), 



(44) 



J — oo 



A{t) = e*^* A{0) e 



-mt 



(45) 



where 7i is the Hamiltonian. Partial integration of (44) yields with the help 
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of (45) the equations of motion 

z{{A;B))i^^ = ([Ai?](±)> + (([AW];S))i±) (46) 
= {[A,B]^^))-{{A;[B,n]))i^\ (47) 

where the ensemble averages () can be related selfconsistently to the respec- 
tive admittance functions (see below). For most of the nontrivial Hamiltoni- 
ans of interest, the equations of motion for the relevant double-time Greens 
functions constitute an infinite hierarchy, whose termination {e.g., in the 
Hartree-Fock approximation) has been a frequently employed approximation 
technique in the theory of magnetism. It is therefore of interest to consider 
models where the exact double-time functions can in fact be evaluated in 
closed form. 

For large frequencies, the double-time functions decrease at least as 
{{A ; B)) ~ z^^, as can be readily seen from their equations of motion (46) 
and (47). Moreover, these functions are analytic off the real frequency (lo) 
axis, across which there exists a cut discontinuity 

{{A ; S)),±,o = {{A ; B))l ± i{{A ; B))l (48) 

where the prime superscript denotes the real part, while the double prime 
stands for the imaginary part of the function. Thus the response functions 
may be represented using the Hilbert transformation 

= (49) 

which is also called spectral representation. 

For the anticommutator (+) metric, the time-correlation functions (ex- 
pectation values) may be obtained using 

{SA{t)SBm = [^e-'-'[l-f{u)]G'U^) (50) 

J IT 

{6B{0)SA{t)) = [^e-^^'f{u^)G'U^), (51) 

where SA = A ~ {A), SB = B - (B), Gab{z) = -{{A;B))+ and /{lo) is 
the Fermi distribution function (temperature units are chosen such that the 
Boltzmann constant is unity). For the commutator (response, (— )) functions 
we have 

{6A{t) dB{0)) = e--* [1 + n(a;)] xabH (52) 

{dB{Q)6Ait)) = [^e-^-'niu)xABi^), (53) 

J TT 
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where the A — B susceptibiUty is defined as xab{z) = —{{A;B))^ and n{uj) 
is the Bose function. Furthermore, the fluctuation-dissipation theorem 



gives a relation between the commutator and anticommutator functions. 



Part of this work was carried out at the Institute fiir Theorie der Kon- 
densierten Materie at the Universitat Karlsruhe, Germany. MMS is grateful 
to P. Wolfle for cordial hospitality, G. Schon for discussions and the HER- 
AEUS Foundation for financial support. MPVS acknowledges the support by 
the Finnish Cultural Foundation. This work has also been supported by the 
Academy of Finland through the project "Theoretical Materials Physics". 



1. R. C. Ashoori, Nature 379, 413 (1996). 

2. S. M. Croncnwett, T. H. Oosterkamp, and L. P. Kouwenhoven, Science 281, 
540 (1998). 

3. N. H. Bondaco, J. Erland, D. Gammon, D. Park, D. S. Katzcr, D. G. Steel, 
Science 282, 1473 (1998). 

4. Y. Meir, N. S. Wingreen, and P. A. Leo, Phys. Rev. Lett. 66, 3048 (1991). 

5. Y. Alhassid, Rev. Mod. Phys. 72, 895 (2000). 

6. P. W. Anderson, Phys. Rev. 124, 41 (1961). 

7. G. Rickayzen, Green's Functions and Condensed Matter (Academic Press, 1981); 
S. Doniach and E. H. Sondheimer, Green's Functions for Solid State Physicists 
(Imperial College Press, 1974); E. N. Economou, Green's Functions in Quantum 
Physics (Springer, 1983). 

8. W. Weller, Phys. Status Solidi (b) 162, 251 (1990). 

9. H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, Phys. Rev. B 21, 1003 
(1980); 21, 1044 (1980). 

10. R. Bulla, A. C. Hewson, and Th. Pruschke, J. Phys.: Condens. Matter 10, 8365 
(1998). 

11. P. B. Wiegmann and A. M. Tsvelick, J. Phys. C: Solid State Phys. 16, 2281 
(1983); 16, 2321 (1983). 

12. P. S. Riseborough, Phys. Rev. B 67, 045102 (2003). 

13. K. Yosida and K. Yamada, Prog. Theor. Phys. Suppl. 46, 244 (1970). 

14. K. Yamada, Prog. Theor. Phys. 53, 970 (1975). 

15. M. M. Salomaa, Solid State Commun. 38, 815 (1981); 39, 1105 (1981). 

16. F. Wei, J.-L. Zhu, and H.-M. Chen, Phys. Rev. B 67, 125410 (2003). 

17. J. Brinckmann, Phys. Rev. B 54, 10465 (1996). 

18. M. E. Foglio and M. S. Figueira, Phys. Rev. B 60, 11361 (1999). 

19. G. Czycholl, A. L. Kuzemsky, and S. Wermbter, Europhys. Lett. 34, 133 (1996). 




(54) 



ACKNOWLEDGMENTS 



REFERENCES 



Response functions of an etrtificial Anderson atom in the atomic limit 



20. D. Meyer, T. Wegner, M. Potthoff, and W. Nolting, Physica B: Condensed 
Matter 270, 225 (1999). 

21. O. Takagi and T. Saso, J. Phys. Soc. Japan 68, 2894 (1999). 

22. A. C. Hewson, The Kondo Problem to Heavy Fermions (Cambridge University 
Press, 1993). 

23. P. W. Anderson, Phys. Rev. 112, 1900 (1958). 

24. The superscript zero denotes static response. 

25. P. W. Anderson, Phys. Rev. Lett. 34, 953 (1975). 

26. S. R. Micnas, J. Ranninger, and S. Robaszkiewicz, Rev. Mod. Phys. 62, 113 
(1990). 

27. A. Taraphder and P. Coleman, Phys. Rev. Lett. 66, 2814 (1991). 

28. R. D. Harris, J. L. Newton, and G. D. Watkins, Phys. Rev. B 36, 1094 (1987). 

29. S. Ghosh and V. Kumar, Phys. Rev. B 46, 7533 (1992). 

30. R. C. Ashoori, H. L. Stormer, J. S. Weiner, L. N. Pfeiffer, S. J. Pearton, K. W. 
Baldwin and K. W. West, Phys. Rev. Lett. 68, 3088 (1992). 

31. R. C. Ashoori, H. L. Stormer, J. S. Weiner, L. N. Pfeiffer, S. J. Pearton, K. W. 
Baldwin, and K. W. West, Physica (Amsterdam) 189B, 117 (1993). 

32. N. B. Zhitenev, R. C. Ashoori, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 
79, 2308 (1997). 

33. R. C. Ashoori, N. B. Zhitenev, L. N. Pfeiffer, and K. W. West, Physica (Ams- 
terdam) 3E, 15 (1998). 

34. M. Brodsky, N. B. Zhitenev, R. C. Ashoori, L. N. Pfeiffer, and K. W. West, 
Phys. Rev. Lett. 85, 2356 (2000). 

35. Y. Wan, G. Ortiz, and P. Phillips, Phys. Rev. Lett. 75, 2879 (1995). 

36. M. E. Raikh, L. 1. Glazman, and L. E. Zhukov, Phys. Rev. Lett. 77, 1354 (1996). 

37. M. M. Salomaa, Z. Phys. B 25, 49 (1976). 

38. J. Hubbard, Proc. R. Soc. A 281, 401 (1964). 

39. J. Hopkinson and P. Coleman, Phys. Rev. B 67, 085110 (2003). 

40. S. C. Bradley, R. Bulla, A. C. Hewson, and G. M. Zhang, Eur. Phys. J. B 11, 
535 (1999). 

41. P. Schlottmann and A. A. Zvyagin, Phys. Rev. B 67, 115113 (2003). 

42. D. N. Zubarev, Sov. Phys. Usp. 3, 320 (1960). 

43. J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 (1966). 



